1 999ASPC . .194. . .44H 


Optical  and  IR  Interferometry  from  Ground  and  Space 
ASP  Conference  Series,  Vol.  194,  1999 
S.C.  Unwin,  and  R.V.  Stachnick,  eds. 


Imaging  and  Modeling  of  Binaries  with  NPOI 


C.A.  Hummel 

US  Naval  Observatory,  3450  Massachusetts  Avenue  NW, 

Washington,  DC  20392,  USA 

Abstract.  This  paper  focuses  first  on  ground  based  observations  with 
long  baseline  optical  interferometers  of  a  class  of  binaries  featuring  double- 
lined  or  variable  spectra.  Following  a  short  perspective  on  implications 
of  this  technique,  we  present  imaging  and  modeling  results  of  two  dou¬ 
ble  stars  (£  Orionis  A  and  o  Leonis)  observed  with  the  Navy  Prototype 
Optical  Interferometer  (NPOI).  Brief  introductions  are  given  on  both  top¬ 
ics,  as  well  as  an  outline  of  current  plans.  Secondly,  we  simulate  future 
imaging  observations  of  stellar  disks  with  NPOI,  and  elaborate  on  related 
issues  of  array  design  and  wavelength  coverage.  We  point  out  differences 
between  optical  interferometry  and  the  well  established  technique  of  radio 
interferometry  on  very  long  baselines. 


1.  Introduction 

In  a  manner  of  speaking,  the  extraordinary  resolving  power  of  modern  optical 
long  baseline  interferometers  is  about  to  lead  to  the  disappearance  of  the  spec¬ 
troscopic  binary  as  a  species.  Not  more  than  15  years  ago,  an  unbridgeable 
gap  existed  between  the  vast  majority  of  stars  exhibiting  periodic  line  shifts  as 
well  as  composite  spectra  and  what  direct  imaging  methods  could  do  to  see  the 
companions. 

These  days,  interferometers  like  COAST  (The  Cambridge  Optical  Aperture 
Synthesis  Telescope,  UK)  and  NPOI  (The  Navy  Prototype  Optical  Interferom¬ 
eter,  AZ)  routinely  produce  images  on  scales  of  a  few  milliarcseconds,  and  are 
about  to  resolve  the  majority  of  the  “spectroscopic”  binaries.  Such  images  have 
been  published  of  Capella  (Baldwin  et  al.  1996)  and  of  Mizar  A  (Benson  et  al. 
1997).  The  latter  star  has  a  long  history  of  astronomical  firsts.  It  was  the  first 
double  star  (components  A  and  B)  to  be  discovered  by  J.  B.  Riccioli  of  Bologna, 
Italy  (and  possibly  even  earlier  by  G.  Galileo,  see  a  web  page  by  L.  Ondra  at 
http://www.bm.cesnet.cz/~ondra/mizar/article.htm),  the  first  to  be  pho¬ 
tographed  by  G.  P.  Bond  of  Harvard  College  Observatory  in  1857,  and  the  first 
spectroscopic  binary  (component  A)  measured  by  E.  Pickering,  also  at  Harvard. 
Now,  a  very  accurate  orbit  of  Mizar  A  has  been  published  (Hummel  et  al.  1998), 
with  median  residuals  of  70  microarcseconds. 

Ground  based  optical  interferometers  feature  long  baselines  of  up  to  several 
hundreds  of  meters  and  achieve  very  high  resolution  despite  atmospheric  turbu¬ 
lence  since  they  are  designed,  in  effect,  as  adaptive  telescopes.  A  good  review 
of  the  current  instrumentation  can  be  found  in  Armstrong  et  al.  (1995). 

44 


©  Astronomical  Society  of  the  Pacific  •  Provided  by  the  NASA  Astrophysics  Data  System 


Report  Documentation  Page 


Form  Approved 
OMB  No.  0704-0188 


Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington 
VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  a  penalty  for  failing  to  comply  with  a  collection  of  information  if  it 
does  not  display  a  currently  valid  OMB  control  number. 


1.  REPORT  DATE 

1999 


2.  REPORT  TYPE 


4.  TITLE  AND  SUBTITLE 

Imaging  And  Modeling  Of  Binaries  With  NPOI 


6.  AUTHOR(S) 


7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

U.S.  Naval  Observatory ,3450  Massachusetts  Avenue, 
N.W., Washington, DC, 20392 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 


3.  DATES  COVERED 

00-00-1999  to  00-00-1999 

5a.  CONTRACT  NUMBER 

5b.  GRANT  NUMBER 

5c.  PROGRAM  ELEMENT  NUMBER 

5d.  PROJECT  NUMBER 

5e.  TASK  NUMBER 

5f.  WORK  UNIT  NUMBER 

8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 


10.  SPONSOR/MONITOR'S  ACRONYM(S) 

11.  SPONSOR/MONITOR'S  REPORT 
NUMBER(S) 


12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  unlimited 

13.  SUPPLEMENTARY  NOTES 

Proceedings  from  ASP  Conference  Vol.  194.  Edited  by  Stephen  Unwin  and  Robert  Stachnik.  ISBN: 
1-58381-020-X  (1999),  p.44 

14.  ABSTRACT 

The  Navy  Prototype  Optical  Interferometer  is  operational  with  three  baselines  of  up  to  38  meters,  which 
can  be  combined  to  form  a  closure  phase  usable  for  obtaining  images.  Work  is  well  underway  to  extend  the 
baseline  lengths  to  a  maximum  of  437  meters,  as  well  as  to  implement  a  6-way  beam  combiner  which  will 
enable  the  simultaneous  operation  of  15  baselines.  Since  NPOI  records  data  over  a  wavelength  range  of  500 
nm  to  850  nm,  the  coverage  obtained  with  Earth-rotation  aperture  synthesis  rivals  the  aperture  coverage  of 
narrow-band  multi-station  radio  interferometers.  In  order  to  investigate  the  technique  of  obtaining  images 
from  NPOI  data,  we  present  actual  imaging  results,  as  well  as  images  from  simulated  NPOI  observations  of 
stellar  surfaces  and  multiple  stellar  systems.  We  discuss  issues  including  the  wavelength  dependence  of 
stellar  structure,  as  well  as  their  variability  during  the  observations.  We  also  point  out  fundamental 
differences  to  the  imaging  of  radio  data. 

15.  SUBJECT  TERMS 


16.  SECURITY  CLASSIFICATION  OF: 

17.  LIMITATION  OF 

18.  NUMBER 

19a.  NAME  OF 

ABSTRACT 

OF  PAGES 

RESPONSIBLE  PERSON 

a.  REPORT 

unclassified 

b.  ABSTRACT 

unclassified 

c.  THIS  PAGE 

unclassified 

Same  as 
Report  (SAR) 

8 

Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std  Z39-18 


1 999ASPC . .194. . .44H 


Imaging  and  Modeling  of  Binaries  with  NPOI  45 

2.  The  Navy  Prototype  Optical  Interferometer 

Conceptually,  the  NPOI  looks  like  the  VLA  (The  Very  Large  Array,  NM,  a  con¬ 
nected  element  radio  interferometer),  because  of  its  Y-configuration  of  stations. 
Each  arm  in  the  Y  is  250  m,  providing  baselines  of  up  to  440  m.  At  optical  wave¬ 
lengths,  this  gives  a  resolution  of  better  than  300  fj,as.  There  are  30  stations  all 
together.  Combined  with  a  vacuum  feed  system,  vacuum  delay  lines,  six  move- 
able  siderostats  featuring  35  cm  clear  aperture,  and  32  channels  covering  450 
nm  to  850  nm  including  IR  capabilities,  the  NPOI  will  be  a  powerful  research 
tool  for  multiple  stellar  systems  and  disks  of  stars.  A  thorough  description  can 
be  found  in  Armstrong  et  al.  (1998). 


3.  Imaging 

The  recent  achievement  of  producing  images  (rather  than  just  models)  of  stars 
has  highlighted  the  potential  of  optical  interferometers.  A  great  deal  of  experi¬ 
ence  dealing  with  imaging  sparse  aperture  interferometric  data  has  been  gathered 
by  VLBI  (Very  Long  Baseline  Interferometry,  intercontinental  baselines  of  radio 
telescopes),  a  radio  technique  now  routinely  producing  high  dynamic  range  maps 
(>  1000)  using  phase  and  amplitude  self-calibration  as  well  as  image  deconvolu¬ 
tion  algorithms.  Whereas  in  VLBI  the  fundamental  observables  are  the  visibility 
amplitude  and  phase,  in  optical  interferometry  they  are  the  squared  amplitude 
and  closure  phase.  Closure  phases  are  the  sum  of  the  visibility  phase  between 
three  (or  more)  stations.  They  are  the  only  phases  free  of  atmospheric  phase 
noise  and  thus  suitable  for  imaging  provided  the  missing  phase  information  is 
substituted  with  model  phases.  This  is  the  essence  of  phase  self-calibration.  In 
addition,  optical  interferometers  can  observe  very  wide  band  passes,  which  can 
be  used  to  fill  in  the  aperture  coverage  in  order  to  enhance  the  stability  of  the 
imaging  process,  but  requires  solving  (or  accounting)  for  the  wavelength  depen¬ 
dence  of  the  source  structure.  Such  imaging  algorithms  have  not  been  developed 
yet. 

For  the  images  in  this  work,  we  have  taken  the  square  root  of  the  squared 
visibility  amplitude,  and  assigned  the  one  closure  phase  (three  station  observa¬ 
tions)  to  one  of  the  baselines  while  setting  the  other  two  baseline  phases  identical 
to  zero.  The  data  was  then  imaged  using  the  AIPS  (Astronomical  Image  Pro¬ 
cessing  System)  software  commonly  used  for  radio  data.  Typically,  the  model 
converged  in  about  five  iterations. 


4.  Modeling 

In  order  to  obtain  astrophysical  parameters  of  the  stellar  systems,  modeling 
is  the  best  and  most  robust  technique  if  the  source  structure  can  be  described 
with  just  a  few  well  defined  structural  elements  (disks,  uniform  or  limb  darkened, 
possibly  with  spots,  etc.).  Data  from  different  techniques  (e.g.,  interferometry, 
spectroscopy,  photometry)  can  be  combined  to  better  constrain  the  model  pa¬ 
rameters.  We  have  developed  a  hierarchical  model  format  including  physical 
parameters,  like  masses  and  effective  temperatures,  for  fitting  the  data.  In  ad- 
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dition,  model  atmosphere  spectra  are  used  to  obtain  fluxes  and  limb  darkening 
coefficients  which  are  integrated  over  the  NPOI  channel  passbands. 

In  the  case  of  binaries,  orbital  elements  are  part  of  the  set  of  model  param¬ 
eters,  and  can  be  fit  either  directly  to  the  visibility  data  if  the  orbital  period  is 
short,  or  to  the  relative  positions  determined  either  from  interferometry  or  by 
other  means  (e.g.,  speckle). 


5.  Close  binaries 

Continuing  a  program  initiated  at  the  Mark  III  interferometer  (e.g.,  Hummel 
et  al.  1995),  observations  of  spectroscopic  binaries  have  been  made  with  NPOI 
in  order  to  resolve  the  components  and  map  their  relative  orbits.  In  addition 
to  measuring  the  relative  component  positions,  the  NPOI  data  also  allow  us 
to  determine  their  magnitude  difference  and  thus  the  individual  components’ 
brightnesses. 


Semi-major  axis  =  4.50  mas 
Inclination  =  57.7  deg 
Line  of  nodes  =  190.3  deg 
Epoch  =  JD  2448223.10 
Period  =  14.498  d 

Orbit  is  circular 


Figure  1.  The  apparent  orbit  and  preliminary  elements  of  o  Leonis 
with  data  from  NPOI  and  the  Mark  III. 


6.  Wide  binaries 

Binaries  with  component  separations  between  several  tens  of  mas  and  about  150 
mas  are  observable  both  with  NPOI  and  speckle  interferometry.  As  the  latter 
method  does  not  provide  good  relative  photometry  of  the  components,  NPOI 
photometry  in  the  V  and  I  bands  adds  important  information  for  deriving  com¬ 
ponent  luminosities.  We  have  selected  from  the  WDS  (The  Washington  Double 
Star  Catalog,  maintained  at  the  US  Naval  Observatory)  systems  with  current 
separations  of  less  than  200  mas,  and  are  observing  some  of  them  routinely.  We 
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present  here  maps  of  the  bright  09.7  supergiant  £  Orionis  A,  in  which  we  found 
a  companion  about  2  magnitudes  fainter  than  the  primary.  The  orbital  motion 
can  clearly  be  seen  and  a  preliminary  orbit  suggests  a  period  of  about  7.5  years. 


PLot  fito  v«r>ton  1  cr*«t*d  O6-APP-1006  16:49:37 
0SC104*  I  POL  3X204E+14HZ  BSC1046.ICLN.33 


c*nl*r  *t  RA  05  40  45.5*714  DEC  -01  56  33.2601 
Pwtkflux-  B.1170E-O1  JY/BEAM 
UV4-  8,1 170E-O3  *  (  -.900,  .500,1.000, 

2.000,  5.000,  10X0,  20.00,  50X0,  00.00) 


PLot  flk  VM-sion  1  erMted  02-MAY-1999 12:12:26 
B  SCI  048  I  POL  4J750E+14HZ  BSC1044JCLN.58 


Center  at  HA  05  40  4&5X714  DEC -01  56  33.2801 

pMk  flux  -  8.O07OE-O1  JY/BEAM 

urn  r  B.097OE-C3  *  (  -.500,  .500,  1.000, 

2X00,  5X00,  10X0,  20X0,  50.00,  8a 00) 


Figure  2.  Maps  of  £  Orionis  A  for  February  1998  and  1999. 


7.  Stellar  masses  and  luminosities 

If  one  selects  from  the  catalog  of  Batten  et  al.  (1989)  all  spectroscopic  binaries 
(with  orbits)  brighter  than  V  =  6,  north  of  —10°,  and  with  estimated  semi- major 
axes  between  1  mas  and  200  mas,  about  300  systems  are  found  with  spectral 
types  covering  the  majority  of  the  HR  diagram.  A  program  aimed  at  system¬ 
atically  observing  such  a  list  of  binaries  with  a  complementary  spectroscopic 
program  to  detect  the  secondaries  could  significantly  increase  the  number  of  ac¬ 
curate  stellar  masses  measured.  The  non-redundant  array  configuration  shown 
below  could  be  used  to  survey  a  large  number  of  systems  for  multiplicity. 


8.  Stellar  surfaces 

The  NPOI  array  was  specifically  designed  to  be  operated  in  different  configu¬ 
rations  in  which  long  baselines  can  be  bootstrapped  by  tracking  the  fringes  on 
shorter  ones  where  the  fringe  contrast  is  high  enough.  This  is  because  in  the 
regime  where  the  longest  baselines  begin  to  see  surface  structure  on  a  stellar 
disk,  the  fringe  contrast  has  dropped  to  very  low  levels  as  the  disk  is  completely 
resolved. 

For  the  imaging  simulation  shown  below,  we  added  Gaussian  noise  of  o  — 
0.05  to  the  real  and  imaginary  part  of  the  simulated  complex  visibilities.  All  32 
channels  of  the  NPOI  were  used.  Two  iterations  of  phase  self-calibration  starting 
with  a  spherical  component  were  performed  while  deconvolving  the  image  with 
a  Maximum  Entropy  algorithm  (AIPS  task  VTESS).  Compared  to  the  widely 
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00:  t  2  3,  01:  1  2  3  4  5  6,  Ch:  1  -32 


U  [MX] 


OB:  1  2  3.  Bl:  12  3  4  5  6,  Ch:  1-32 


U  [MX] 


Figure  3.  Earth- rotation  aperture  synthesis  with  a  non-redundant 
survey  configuration  (left)  and  a  redundant  stellar  disk  imaging  con¬ 
figuration  (right). 


OB:  1  2  3.  Bl:  12  3  4  5  6,  Ch:  1  -32 


UV-RodiuS  [MX] 


00:  1  2  3,  Bl:  1  2  3  4  5  6.  Ch:  1  -32 


Figure  4.  The  squared  visibility  amplitude  as  a  function  of  uv-radius 
for  a  12  mas  stellar  disk  observed  with  the  redundant  array  configura¬ 
tion.  Note  the  very  low  amplitudes  at  spacings  with  the  most  informa¬ 
tion  of  surface  structures  on  the  disk. 
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Figure  5.  Image  of  a  stellar  spot  using  simulated  data  from  NPOI 
in  the  redundant  array  configuration.  The  stellar  disk  is  12  mas  in 
diameter,  the  spot  diameter  is  3  mas,  the  restoring  beam  is  2  mas. 
The  spot  has  a  temperature  of  3500  K,  whereas  the  star  is  a  5000  K 
limb  darkened  disk. 


used  CLEAN  algorithm,  the  former  is  better  suited  for  the  large  and  low  contrast 
area  represented  by  the  stellar  disk. 
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Discussion 

Andrew  Gould:  How  long  until  you  measure  those  164  stars  with  1%  mass 
precision? 
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Christian  Hummel:  Based  on  the  Batten  et  al.  Catalog  of  spectroscopic  binary 
elements  I  selected  all  stars  brighter  than  V=6,  north  of  -10  degrees  declination, 
and  with  estimated  semimajor  axes  larger  than  1  mas,  and  derived  164  for  the 
number  of  components  whose  masses  could  be  determined  with  better  than  1% 
precision  if  the  relative  astrometry  had  a  precision  of  10  micro-as  and  the  velocity 
amplitudes  a  precision  of  50  m/s.  These  are  difficult  to  achieve,  but  maybe  not 
impossible.  Such  a  project  would  take  many  years. 

Vassilis  Zarifis:  What  is  the  effect  of  the  missing  points  in  the  (u,v)  plane  on 
the  synthetic  imaging? 

Christian  Hummel:  Missing  UV  -  points  cause  sidelobes  around  bright  sources 
in  the  image.  They  can  be  effectively  removed  with  the  CLEAN  algorithm. 
Large  gaps  in  the  UV  plane,  however,  will  degrade  the  image  so  much  that  only 
the  simplest  of  source  structure  can  be  imaged. 


Slava  Turyshev,  Stephen  Unwin  and  Rudolf  Danner 
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